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The response function of domain growth processes, and in particular the violation of the fluctuation- 
dissipation theorem, are studied both analytically and numerically. In the asymptotic limit of large 
times, the fluctuation-dissipation ratio X, which quantifies this violation, tends to one if C > m 2 
and to zero if C < m 2 , corresponding to the fast ('bulk') and slow ('domain- wall') responses, 
respectively. In this paper, we focus on the pre-asymptotic behavior of the domain-wall response. 
This response is shown to scale with the typical domain length L(t) as 1/L(t) for dimension d > 2, 
and as ln(L(t)) / L(t) for d — 2. Numerical results confirming this analysis are presented. 
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I. INTRODUCTION 

Domain growth systems are the paradigm of systems 
that do not reach equilibrium. Hence, it has been a re- 
curring theme in the field of spin and structural glasses 
to think these systems as displaying some form of coars- 
ening . In such non-equilibrium systems, time trans- 
lation invariance does not hold, and all time dependent 
correlation functions and response functions depend on 
two times (the time origin corresponding generally to the 
time at which the system has been quenched into the non- 
equilibrium state). In domain growth phenomena, an au- 
tocorrelation function C(t,t') of the form C(L(t)/ L(t')) 
is usually interpreted as arising from structures whose 
size grows as L(t). A similar functional form, however, 
is also found to describe the out of equilibrium dynamics 
of mean-field models of glasses || , although there is by 
construction no length L in such models. In fact, the 
difference between both kinds of models only becomes 
manifest when one also considers the response functions 
associated to the correlation functions. 

Generally speaking, the fluctuation-dissipation theo- 
rem (FDT), which for equilibrium systems relates the 
response functions to the correlation functions, does not 
hold in systems that are out of equilibrium. The de- 
viations from FDT are conveniently described by intro- 
ducing the Fluctuation Dissipation Ratio X(t, t') defined 
through TR(t,f) = X(t,t')^(t,t'), where R is a re- 
sponse function and C the associated correlation func- 
tion. In mean field models of spin glasses, the behav- 
ior of this FDR has been well established, at least in 
the asymptotic limit of large times. One has X ~ 1 for 
"fast" processes (t ~ t', large C), but a value < X < 1 
for the "slow" processes corresponding to well-separated 
times @J3| . This observation prompted several work- 
ers [!U^] to calculate the large-time response of pure (non- 



disordered) coarsening models, in order to quantify the 
similarities and differences with mean field models of 
glasses. In order to make a comparison with the glassy 
case, one computes the staggered response to a spatially 
random field, to make the perturbation uncorrelated with 
the equilibrium pure states (as is the case, for example 
of a uniform field for a spin glass) . 

The result is that X ~ for all but the smallest time 
differences. In other words, the long-time memory of 
coarsening systems tends to vanish, unlike in mean-field 
glass models where it does not since X > even at long 
times. As far as we know, all systems in which two (or 
only few) phases separate have X = at long times, 
and this has been proven under certain assumptions [Q. 
Physically, this feature can be understood from the fact 
that for long times the response will be dominated by 
the bulk response of the domains that form during the 
coarsening process. The response at time t of a spin to 
a field applied at time t' will be nonzero only if the spin 
is not swept by a domain wall between t and t' . Other 
types of response involve the domain walls themselves, 
whose density decreases with time, and therefore vanish 
in the limit of large times. 

From the experimental point of view, aging experi- 
ments H show that glasses such as spin glasses or molec- 
ular glasses do have long term memory. The asymptotic 
nature of experimental results is however, always ques- 
tionable. This is even more the case in numerical stud- 
ies, which for a number of models (spin glasses, struc- 
tural glasses, kinetic models and polymers in disordered 
media |J) have obtained results in qualitative agree- 
ment with mean-field theory (X ^ 0). It is therefore 
a relevant question to study the deviations from FDT 
in the pre-asymptotic limit. An understanding of this 
pre-asymptotic behavior should allow to distinguish be- 
tween true long term memory and a slow approach to a 
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vanishing X. 

In this paper we present such a study for the ferro- 
magnetic coarsening or phase-separation of pure (non- 
disordered) systems after a quench at time t = from 
a homogeneous phase (T = oo) into a two-phase re- 
gion (T < T c ), with and without local conservation of 
the order parameter. In the thermodynamic limit, the 
equilibrium state, where the two phases are completely 
separated is never achieved. We confirm the previous re- 
sults for the absence of long term memory in the response 
function , and then study the scaling of this response 
in the pre- asymptotic regime (large but finite times). It 
turns out that the model-dependence enters only through 
the form of the growth law L(t). 

We present the systems and the dynamical quantities 
under study in Section O. The numerical simulation is 
described in the Section jnj, while the analytic study is 
presented in the Section IV. 



II. MODELS AND DEFINITIONS 

The systems considered here will be described by a 
coarse-grained formulation, with a scalar order parame- 
ter </>(r, t) and a Ginzburg-Landau free energy functional 
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where h(r,t) is the field conjugated to <f>(r,t). Experi- 
mental situations under consideration are for example the 
coarsening in a ferromagnet, or spinodal decomposition 
in a binary alloy. 

Domain growth processes have been much studied 
since the early works of Lifshitz, Slyosov and Wagner. 
Ref. [jlO| is a very complete review on the topic. 

If the order parameter is not locally conserved, we have 
the Ginzburg-Landau equation 
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dt 



(2.2) 



where 77 is a gaussian markovian noise term satisfying 
(r](r, t)) = and (rj(r, t)rj(r', if)) = 2TS(r - r')5(t - if), 
and T is the temperature. When the order parameter 
is conserved, the evolution is given by the Cahn-Hilliard 
equation 



dt 



SF\ 



(2.3) 



In that case, the thermal noise is characterized by the 
two moments of the gaussian distribution (r](r,t)) = 0, 
and (r)(r, t)r](r', if)) = -2T6(t - t')V 2 5(r -r'). 

Interesting dynamical quantities in the study of the out 
of equilibrium properties are the autocorrelation function 
defined by 

C(t, t w ) = ^J d d r (<t>(r, t)<j>{r, t w )) , (2.4) 



and the associated response function R(t,t w ) = 
(6(j)(t)/Sh(tw)) . At equilibrium, these two quantities de- 
pend on time difference r = t — t w only, and are related 
by the usual fluctuation dissipation theorem 



R(r) 
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T dr ' 

Out of equilibrium we write jjj: 

X(t, t w )dC{t,U 



R(t, t w ) — 



T 



(2.5) 



(2.6) 



which defines X(t, t w ) as the fluctuation dissipation ratio. 

The strategy for our study is now standard. The 
quench of the system takes place at t = 0. In order 
to compute the correlation C(t,t w ) we record the con- 
figurations of the system evolving at zero external field, 
h = for times t > t w . The operation is repeated on 
several samples in order to improve the statistics. 

The integral of the linear response function M(t, t w ) = 
f t dsR(t,s) is computed by letting the system evolve 
under the influence of a small field switched on at t w , 
and recording the magnetization at time t. The field is 
random in space and stationary [||. It is drawn from 
a gaussian distribution with first moment h(r) — and 
second h(r)h(r') — ho 2 S(r — r'), respectively. In the 
language of magnetic systems, the integrated response 
function is thus the staggered magnetization 



M(t,t w ) 



1 



h 2 V 



d d r(h(r)4>(r,t)). (2.7) 



An important property of the FDR has to be empha- 
sized for the present discussion. This property, found 
analytically within mean-field models and verified nu- 
merically in various glassy systems ^^], is that in the 
asymptotic regime of t, t w — > 00, X(t, t w ) depends on the 
times only through a non singular function of correlation 
function C(t,t w ), that is X(t,t w ) = x{C(t,t w )). When 
this property holds, the generalized FDT ( |2.6| ) gives the 
following relation between M(t, t w ) and C(t, t w ): 



I r C{t,t) 
M(t,t w ) = - / dCx{C). 



T 



C(t,t w ) 



(2.8) 



In equilibrium systems, x = 1, so that one has the 
relation TM(t,t w ) = C(t,t) - C(t,t w ). More gen- 
erally, in non-equilibrium systems, a parametric plot 
M(t,t w ) vs C(t,t) — C{t,t w ) is independent of t w , and 
allows a direct determination of x(C). 

In the pre-asymptotic regime, the parametric plot of 
M(t,t w ) vs C(t,t w ) (with t as the parameter) will gen- 
erally depend on t w . Interesting information can never- 
theless be extracted from this plot, as will be seen in the 
next section. In particular, a constant slope is indicative 
of a constant value of X, and a zero slope (plateau in M) 
corresponds to a loss of memory in the response. 
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A second property of the FDR is that under certain 
assumptions H it happens to coincide with the static 
Parisi function x(q) = dq' P(q'), where P(q) is is the 
probability distribution of overlaps between real replicas 
of the same system. For the ferromagnetic case, P(q) is 
trivial and P(q) = 5(q - M 2 ), where M = M(T) is the 
magnetization. Therefore, we expect the FDR to be 1 if 
1 > C > M 2 , and if M 2 > C. 



III. SIMULATION OF A SPINODAL 
DECOMPOSITION 

The Monte-Carlo studies of Ref. |6| agreed qualita- 
tively with the above behavior of the FDR, but to our 
knowledge no quantitative results are available yet. It 
is moreover clear that the asymptotic regime where the 
parametric plot Response/Correlation is supposed to col- 
lapse into a master-curve was not reached, the plot still 
conserving a dependence on t w . In this work, we will be 
interested in a quantitative study of this pre-asymptotic 
behavior. 

For thi s p urpose, the stochastic partial differential 
equation ( |2.3| ) was numerically solved, in order to model 
a spinodal decomposition. Both time and space were 
discretized. A 1024 x 1024 square lattice with periodic 
boundary conditions was used. Spatial derivatives are 
treated using an implicit spectral method. Time deriva- 
tives were approximated using a simple Euler scheme. 
No real improvements have been obtained using a sec- 
ond or fourth order stochastic Runge-Kutta algorithms. 
The following recurrence relation is then obtained in the 
discretized Fourier space: 
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^ tn+l)= [l + {k±-k 2 )At 
x \<f>(k,t n ) - k 2 M(4> 3 (k,t n ) - h(k)) + Vkrj 



(3.1) 



After discretization, the noise term 77 is characterized by 
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The algo- 
rithm is the following. Knowing the fields (j>(r,t n ) and 
<f> 3 (r,t n ), the Fourier transforms 4>(k,t n ) and 3 (fc,t„) 
are computed. The recurrence relation (3.1) is then used 
to obtain <fi(k,t n +i). Fourier transforming again gives 

The influence of the parameters Ax, Ay, T and At 
on the numerical integration is discussed in the litera- 
ture pTp . We chose Ax = Ay = 0.5, in order to get 
mesh-size independent results. The thickness of the do- 
main walls in the late stage of the phase separation is 
indeed about £ = 1/V2, where £ is the correlation length 
of the model ( |2.1| ) . Their structure is hence sufficiently 
well described by the above discretization. 

The role of At is made less crucial by our choice of an 
implicit algorithm. The linear stability analysis of our 
algorithm gives indeed the following results: the "tan- 
gential bifurcation" |ll[ , that is the small k instability, is 



still present, but it is of course physically essential. On 
the contrary, the "subharmonic bifurcation" [jllj does not 
exist any more. Hence, the only restriction on the time 
step is the average magnitude of the noise which has to be 
kept small in order to avoid numerical divergences. We 
chose then the highest possible value of At compatible 
with the temperature T. A small temperature allows a 
large At, but obliges to work with a very small magnetic 
field when computing the response function (see below). 
We chose finally T = 0.1 and At = 0.2 in order to ex- 
plore a large time range. For very short times, a spurious 
behavior (see the caption of the figure ||) related to this 
rather large value of At can be observed. We checked 
however that this deviations vanish when At is smaller, 
and do not affect the long time evolution we are inter- 
ested in. 





FIG. 1. Field configurations during the simulated phase 
separation for times 317, 1262, 5024, 20100. Each color rep- 
resents one phase. 

Fig.|l| presents different field configurations during the 
coarsening process. Looking at these pictures, it is clear 
that the coarsening process can be characterized by the 
typical size L(t) of the domains. Growth laws are well 
known jl0| and are L(t) ~ t 1 / 3 in the conserved case, 
and L{£) ~ t 1 ! 2 in the non-conserved case. The pre- 
ceding remark has a very interesting consequence which 
is known as the scaling hypothesis. As L(t) is the only 
physically relevant length scale, statistical properties of 
the system are the same if we scale all the lengths by the 
factor L(t) @. 

As in experiments, the measure of the domain size is 
obtained by computing the structure factor 



S(k,t) = {cf>{k,t)(j>{-k,t)). 



(3.2) 



The scaling hypothesis implies that it can be written as 
S(k,t) — L(t) d g(kL(t)), where g is a scaling function. A 
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convenient way of obtaining the growth law is to perform 
a circular average of S(k, t) and to compute then (k) = 
f dkS(k,t)k/ J dkS(k,t), which scales as 1/L(t). As in 
Rcf. we obtained the growth law L(t) ~ t 1 / 3 , which 
is valid after a short transient period. The time evolution 
of the circularly averaged structure factor is depicted in 
the inset of the figure @. 
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FIG. 2. Inset: circularly averaged structure factor for 15 
different times from 317 to 23318. The largest time corre- 
sponds to the highest maximum. Main picture: the 15 curves 
collapse on a single scaling function. 

It has a clear maximum, corresponding to the wave- 
vector 27T /L(t). This maximum shifts towards the small 
k, while its amplitude grows with time. The scaling hy- 
pothesis is verified plotting S(k, t)/L(t) 2 vs kL(t), all the 
curves collapsing on a very well defined scaling function 
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FIG. 3. Correlation function for waiting times (from left to 
right) t w = 317, 502, 796, 1262, 2000, 3170, 5024, 7962 and 
12619. Inset: the same curves as a function of t/t w . 

The autocorrelation functions for various waiting times 



t w are shown in the figure [| More precisely, the quan- 
tity C(t, t) — C(t, t w ) is computed, in order to avoid nor- 
malization problems of the correlation functions. As we 
are working with a soft spin system rather than a more 
usual Ising model, C(t, t) is indeed a slowly varying func- 
tion. More precisely C(t,t) ~ C eq — L /L(t), where C eq 
is the equilibrium correlation function in a bulk system 
and Lo/L(t), is proportional to the interface density. 

Two regimes may be distinguished: for times t <C t™, 
the correlation is time translation invariant (TTI), and 
for times t > t w , aging is evident, with the TTI break- 
down, and the correlation falls to 0. This scenario has 
been called weak ergodicity breaking The fluctuation 
dissipation theorem holds in the former, but is violated 
in the latter regime of times. 

It is useful to use this behavior to introduce 



C(t,t w ) — C s t(t t w ) + C ag (t,t w ), 



(3.3) 



where C s t and C ag describe respectively a stationary and 
an aging part in the correlation. 

The scaling hypothesis may be used to predict a scaling 
form for the aging part of the correlation function [id] : 



f 



L(t) 
L(t w ) 



f 



(3.4) 



/ being a scaling function. Eq.(3.4) retains an explicit 
dependence on both times t and t w , typical of an ag- 
ing system. Such a scaling in the correlation function is 
called simple aging. As shown in the inset of figure |^, 
this scaling form describes our results extremely well. 

In order to complete the study of the fluctuation dissi- 
pation theorem, we have to compute the response func- 
tion M(t,t w ) to the static random field applied between 
t w and t (recall Section |ll]) . The field amplitude has to 
be small to obtain a linear response. The best numeri- 
cal test we found for this purpose is the comparison of 
the time evolution of (k) with and without the magnetic 
field. When the field is present, the domain walls may be 
slowed down and even pinned if the field is too strong, 
so that the coarsening process is perturbed. This test 
is very sensitive, and we worked with small field ampli- 
tudes (between ho = 0.035 and 0.09) to ensure that the 
coarsening process was not affected. 

With the correlation and the response functions, the 
parametric plot of M(t, t w ) vs C(t, t w ) may be built. The 
data are shown in figure |], for various values of t w , t being 
the parameter. The curves are averaged over 14 to 32 
realizations of the magnetic field. They are qualitative! 
the same as in the previous Monte-Carlo simulations 
with a first part in which the FDT holds, corresponding 
to the times t <^t w . In a second part, which corresponds 
to times t > t w , the FDT is obviously violated, with M 
having a quasi-horizontal plateau. As discussed above, 
this plateau indicates the loss of long term memory in the 
response (X = at long times), consistent with previous 
expectations. 
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exactly as we did for the correlation function. 
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FIG. 4. Test of the fluctuation dissipation theorem. Wait- 
ing times are t w = 317, 502, 796, 2000 and 5024. The horizon- 
tal line is the equilibrium value of the response and the dashed 
line is the FDT regime with slope 1/T. This line has an in- 
tercept which is slightly positive: this is a time discretization 
effect, which can be made to vanish by reducing At. 

The equilibrium value of the response function has 
been numerically computed by performing a similar simu- 
lation in an homogeneous system. This value is indicated 
in figure |] by an horizontal line (Note that this line can 
also be determined analytically, as discussed in the next 
section). From the arguments presented in the intro- 
duction, it could be expected that the long time plateau 
of the integrated response corresponds to that of a sin- 
gle domain. It is clearly seen from the data in figure ||, 
however, that the approach to this asymptotic value is 
extremely slow. 

As the extra response (i.e. the difference between the 
plateau value of M(t,t w ) and the bulk response) can 
be attributed to the domain wall response, it is tempt- 
ing to try to relate this response to the domain wall 
size. Between waiting times t w = 317 and t w — 5024, 
the size of the domains increases multiplied by a factor 
(5024/317) 1 / 3 - 2.51, while the extra response is divided 
by only 1.36. The pre-asymptotic behavior of the FDR 
seems then to be related to L(t w ) only through a non- 
trivial relation, which we explicitly discuss in the follow- 
ing section. 



IV. ANALYTICAL STUDY OF THE 
FLUCTUATION DISSIPATION RATIO 

In order to study the fluctuation dissipation theorem, 
we have to compute sepa rate ly the correlation and the 
response functions. Eq.(3.4) will be sufficient for the 
present discussion, as we only need a scaling form for 
these functions. The response function may also be split 
into 



We compute first M eq = lim t 



M eq (t—t w ), which 



is in fact the static equilibrium response function of a 
single domain. It may be evaluated exactly at T = 
(within the gaussian approximation) and corrected per- 
turbativcly in powers of T. One easily finds for T = 



M e 



d d k 



1 



fc 2 + l/£ 2 



(4.2) 



Recall that we have £ = l/y/2. In the simulation, the 
space is discretized, and this integral becomes then a dis- 
crete sum over the first Brillouin zone. A numerical eval- 
uation of this sum yields a result in perfect agreement 
with the simulation result obtained for an homogeneous 
system, as described in the previous sect ion. The first 
temperature corrections to equation (4.2) can be com- 
puted exactly (see figure |5|), and are indeed found to be 
negligible at the temperatures we used. 





FIG. 5. The three Feynman graphs representing the "one 
loop^ (proportional to the temperature T) corrections to 



Eq.(4.2) 



M(t,t w ) = M eq {t - t w ) + M ag {t,t w ), 



(4.1) 



Note that the integral (4.2) is divergent in the contin- 
uous theory: one has to introduce a spatial cutoff a, sim- 
ulating the underlying lattice spacing. A convenient way 
of doing this is to multiply the integrand by exp(— k 2 a 2 ). 
Thus, the equilibrium response function scales with the 
cutoff as M eq ~ a 2 ~ d for d > 2, and as M eq ~ ln(a/£) for 
d = 2. 

Next, we compute the "aging part" of the response 
function, which involves the response of the domain walls. 
This can be done using one of the "approximate theories 
for scaling functions" [M , which attempt to give an an- 
alytical expression for the scaling function g of the struc- 
ture factor, or equivalently for its Fourier transform. The 
spirit of these theories is to replace the field cf>(r, t), which 
at the late times of the coarsening process is ±1 outside 
the domain walls, by an auxiliary field m(r,t) varying 
smoothly in space. This allows to derive evolution equa- 
tions for m that may -with further approximations- be- 
come tractable. This method has al read y been used to 
justify analytically the scaling form (3.4) of the correla- 
tion function jl2) , and will be used below for the response 
function. 

Unfortunately, such schemes have not been successfully 
developed for the conserved case. Our strategy will be 
then to study the non-conserved case within these ap- 
proximations, and to give physical arguments to extend 
the validity of our result to the conserved case. For this 
purpose, we generalize a calculation of Bray in reference 
|0 of the response to a uniform magnetic field, to the 
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staggered response to a random field. Following Bray, 
the equation for the auxiliary field m reads (at T = 0) 

— — = V 2 77i — 7i a 7i(,V a Vbm + h\Vm\. (4.3) 
at 

The field h depends now on space. Further simplifica- 



tions, the validity of which we do not discuss here [ 10|,|f3| , 
are the replacements n a nb —* 5 a b/d (circular average), 
and |Vm| — > ((Vrn) 2 ) 1 / 2 . This computation scheme 
is near in spirit of the theory of Ohta, Jasnow and 
Kawasaki |L4|. With these two assumptions, (43) be- 
comes 



dm , , 

— = DV 2 m + a(t)h, 



(4.4) 



with D = (d—l)/d and a(t) = ((Vm) 2 ) 1 / 2 . The solution 
for m is: 



m(k, t) = m(k, 0)e 



-k A Dt 



h{k) J dt'a(t')e 



-k 2 D(t~t') 

(4.5) 



Random initial conditions are conveniently chosen from 
a gaussian distribution with mean zero and variance 
A. The quantity we want to compute is the staggered 
magnetization M ag (t,t w ) = (h(r)(f>(r,t)}/ho 2 . At late 
times of the coarsening process, the usual approximation 
cj> ~ sgn(m) can be made. Using the fact that the fields 
m(k,0) and h(k) are gaussian, it is straightforward to 
obtain the staggered magnetization in term of the field 
m and one gets: 



M ag {t,t v 



2 (h(r)m{r)} 



(4.6) 



Using Eq.(0), the relations (m 2 ) ~ A/(Dt) d/2 and 
a(t) ~ vA / 't( d + 2 )/ ' 4 are also obtained and finally: 



M ag (t,t w ) ~ / d d kM(k,t,t w ) 



M(k,t,t w )= / dt'e-* 2 ^*-*') 



{Dt) d / A 
(£><') (d+2)/4 ' 



(4.7) 
(4.8) 



This integral over k is divergent for large k. As in the 
equilibrium case we introduce a cutoff length a via a term 
exp(— k 2 a 2 ). The integrals can now be performed and 
yield 



1 f(± 



M ag (t, t w ) 

_ ln(tL /2 /a) f / t 
~ ln(a/C)d /2 \*» 
The scaling function _F is given by 



F(A) = lim 



d > 2, 

d = 2. (4.9) 
(4.10) 



a-,o J 1 A 2 -* A'( d + 2 )/ 4 (A-A' + A 2 )^/ 2 



Except in dimension d — 2, the cutoff a disappears if 
the non-equilibrium response is measured in terms of the 
equilibrium one. 

The meaning of this result can be better understood by 
considering the response associated to each spatial length 
scale separately. Defining M(k,t,t w ) as the response to 
a sinusoidal perturbation with wave- vector k, we can dis- 
tinguish between two cases. For a wavelength larger than 
the domain size, k <C 1/L(t w ), we obtain 



M(k,t,t w ) - t d/4 



dt' 



t ,(d+2)/i 



L{t w )G 



t 



(4.11) 



This is precisely Eq.(114) in Ref. |1^]. The response of 
long wavelength modes grows with time, although their 
effect becomes negligible because the number of modes 
with k <C 1/L(t w ) decreases with time. On the other 
hand, for short wavelengths, k ^> 1/L(t w ), the integral 
can be approximated to find 



M(k,t,t w ) 



1 

fc 2 " 



1 

7T/2' 



(4.12) 



which is in fact a very simple result. The susceptibility of 
an elastic surface (a flat domain wall) when a field with 
wave- vector k is applied is proportional to 1/fc 2 , and the 
density of interfaces is proportional to 1/L(t) <~ 1/t 1 / 2 . 

The behavior of the aging part of the response function 
may now be simply evaluated as the sum of two terms. 
The first one is the contribution of small wave-vectors 
k <C 1/L(t w ). We have already shown that for the non- 
conserved case, it became negligible as L — > oo. We 
can safely assume that this is a general statement, as 
the influence of the long wavelengths is even smaller in 
the conserved case (c/the Cahn-Hilliard equation). The 
second one, corresponding to wave- vectors k ^> 1/L(t w ) 

scales as f*f? d d k/k 2 L. The long-time response 



Mag(t, t w ) 



M, 



-F 



c.q 



Lit) 

L(t w y \L(t w ) 

\n(L{t w )/a 
ln(a/QL(t u 



F 



Lit) 
L(t w ) 



d>2, 

d=2, (4.13) 



is entirely dominated by the short wavelengt hs. For a 
non-conserved dynamics, with L ~ i 1 / 2 , Eq.(4.9) is re- 
trieved. 

The result (4.13) is now in a form independent of 



the dynamics (conserved or non-conserved order param- 
eter) of the sy stem . Hence, for any coarsening system, 
Eqs.(3T) and ( 4.13 ) give an analytical evaluation of the 
response associated to the domain walls ratio in the aging 
part. As expected, this response vanishes at long times, 
so that in the asymptotic regime t,t w — > oo, the value 
x(C) — is obtained in all dimensions. 
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O) 

to 



7 

\nM). 



FIG. 6. Test of the scaling (1.13) which predicts a linear 



dependence of L(t w )M ag w.r.t. ln(t w ). The dashed line fits 
this dependence very well. 



The scaling ( 4.1 3| ) is tested in the figure |6[ where we 



plot t^^Mag vs ]n(t w ). Here M ag is defined as the dif- 
ference between the plateau value obtained in the simu- 
lation and the equilibrium response. Our data are obvi- 
ously consistent with the above assumption. Hence we 
conclude that the extra response obtained in the simu- 
lation actually corresponds to the domain wall response, 
and will asymptotically vanish. This vanishing, however, 
is extremely slow, so that we can hardly expect to see it 
in any numerical simulation. 



V. CONCLUSION 

As mentioned above, the important measurable (dy- 
namical) difference between mean-field glass models and 
coarsening models is the presence -or absence- of a long 
term memory in the response functions. In view of the 
relations this has in certain cases with the Parisi function, 
the question as to whether real glasses have a value of X 
that stays different from, or tends slowly to zero is some- 
times taken as the modern version of the old "droplet" 
versus "mean- field" debate of the 80 's. Experimentally, 
however, the difference between X tending logarithmi- 
cally to zero or staying constant might not be very dra- 
matic. What is important, however, is that if one has a 
law for the integrated response of the form 



M ag (t,t w ) ~ A(t w )F 



L(t) 
L{t w ) 



(5.1) 



then if A(t) ~ (L(£)) -1 a large aging response is neces- 
sarily linked to very slow scaling laws L(t). The fact that 
we have found for d — 2 & relation A(t) ~ \n(L(t)) / L(t) 
shows that indeed it is possible for the response to fall 
slower than the inverse of the rate of growth L, and one 
can have relatively large long term memories together 
with rather fast growth laws. 



Another important conclusion of the present study, 
which confirms earlier numerical work, is that the domain 
walls can have a large contribution to the response in the 
pre-asymptotic regime, but almost exclusively given by 
their deformation on relatively short lengths. This elastic 
contribution can be considered thermalized, and its con- 
tribution makes longer the segment of slope 1/T in the 
parametric plot of the integrated response versus correla- 
tion. Apart from that, the plot is flat for smaller values 
of C. This is still very different from mean-field glass 
models, in which the out of equilibrium contribution is 
'thermalized' at an effective temperature different from 
T. 

This work was supported by the Pole Scientifique de 
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